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ABSTRACT 

Atmospheric mass loss from Hot Jupiters can be large due to the close proximity of these planets to their host star 
and the strong radiation the planetary atmosphere receives. On Earth, a major contribution to the acceleration of 
atmospheric ions comes from the vertical separation of ions and electrons, and the generation of the ambipolar 
electric field. This process, known as the “polar wind,” is responsible for the transport of ionospheric constituents 
to Earth’s magnetosphere, where they are well observed. The polar wind can also be enhanced by a relatively 
small fraction of super-thermal electrons (photoelectrons) generated by photoionization. We formulate a simplified 
calculation of the effect of the ambipolar electric field and the photoelectrons on the ion scale height in a generalized 
manner. We find that the ion scale height can be increased by a factor of 2-15 due to the polar wind effects. We 
also estimate a lower limit of an order of magnitude increase of the ion density and the atmospheric mass-loss rate 
when polar wind effects are included. 

Key word: planets and satellites: atmospheres 


1. INTRODUCTION 

In the past two decades and in particular following the 
Kepler mission, hundreds of exoplanets have been detected (e.g., 
Schneider 1995; Mayor et al. 2003). Many of these planets 
are gas giants observed at an extremely close orbit of less 
than 0.1 AU from their host star (an orbital period of less 
than 10 days), and are classified under the term “Hot lupiters” 
(HJs). The unexpected close-in orbit of HJs has stimulated many 
science investigations regarding their formation, evolution, and 
tidal interaction (e.g., Papaloizou et al. 2007 and references 
therein), their magnetic interaction with the host star (e.g., 
Cohen et al. 2010 and references therein), and the structure 
and dynamics of their atmospheres (e.g.. Showman et al. 2008 
and references therein). 

In such a close orbit (especially if the star and the planet 
are tidally locked), His are expected to receive extremely large 
amounts of stellar X-ray and EUV radiation (Penz et al. 2008; 
Cecchi-Pestellini et al. 2009). It has been argued that this high 
EUV radiation can lead to strong photo-evaporation of the plan- 
etary atmosphere and high mass-loss rates (Lammer et al. 2003; 
Baraffe et al. 2004, 2006), leading to less massive planets. How- 
ever, this could not be supported by the observed mass distri- 
bution (Hubbard et al. 2007). Observations of Lya emission 
from the HD209458 system have suggested that the planet oc- 
cupies an inflated hydrogen corona with outflow velocities of 
50-100 km s 1 and a mass-loss rate of about 10 u) g s _1 . How- 
ever, there is a debate on whether the observations are effected 
by the host star or whether the observed features are of planetary 
origin (Vidal-Madjar et al. 2003; Ben-Jaffel 2007; Vidal-Madjar 
et al. 2008). A more recent observation of the system, as well 
as of the HD 189733 system, revealed a smaller mass-loss rate 
of about 10 8 g s' 1 (Lecavelier Des Etangs et al. 2010; Linsky 
etal. 2010). 

On the theoretical side, several models for atmospheric 
escape from His have been developed in recent years. Detailed 
models for the chemistry, photoionization, and aeronomy of HJs 
were developed by Yelle (2004) and by Garcia Munoz (2007). 
Tian et al. (2005) and Murray-Clay et al. (2009) performed 


hydrodynamic calculations of thermally driven atmospheric 
escape, and Stone & Proga (2009), Trammell et al. (201 1), and 
Adams (20 1 1 ) included the planetary magnetic field geometry, 
which confines the escaping gas to regions of open field lines. 
The models above predict mass-loss rates not higher than 
10 I() g s' 1 . Some of the models also included the incoming 
stellar wind and found that the planetary outflow is ought to be 
suppressed by the wind. None of the models predict a sufficiently 
high mass-loss rate so that the planet can be evaporated in a 
relatively short timescale. 

In Earth’s upper atmosphere (as well as in other planets), there 
is a well-observed physical process which plays an important 
role in the acceleration of ions. The polar wind (Banks & Holzer 
1968) is the outflow of planetary ions along open field lines. The 
main driver for this process is the ambipolar electric field, which 
is proportional to the electron pressure gradient. Since electrons 
are more mobile than ions, a charge separation is created along 
the magnetic field direction, leading to an electric potential that 
acts on the ions to retain charge neutrality. The end result is an 
acceleration of the ions by this electric field so that the ions are 
dragged by the electrons. 

The electric field applies a force proportional to the negative 
gradient of the electron pressure. Using some simplifications, 
the resulting force is approximately equivalent to half of the 
gravitational force on the major ion species and directed op- 
positely. Since 0 + is the major ion species in Earth’s upper 
ionosphere, the result is a supersonic flow of H + and an increase 
in the 0 + scale height. In addition, photoelectrons, which are 
electrons highly energized due to photoionization (the tail of 
the distribution function), can significantly increase the electron 
temperature, leading to an enhancement of the ions acceleration 
(Lemaire 1972). In Earth’s upper atmosphere, the velocity of 0 + 
is lower than the escape velocity. Nevertheless, 0 + is observed 
to serve as a significant plasma source in the magnetosphere 
(Lennartsson et al. 1981). Tam et al. (1995, 1998) have demon- 
strated by numerical simulation that photoelectrons indeed can 
accelerate 0 + and H + , though they obtained an unrealistic elec- 
tron temperature of 40, 000 K. An additional simulation by 
Khazanov et al. (1997) resulted in a more realistic electron 
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temperature of 16,000 K. Recent numerical simulations by 
Glocer et al. (2012) also included the effects of photoelectrons 
to look at the global outflow solution and compared with in 
situ observations. Their simulations showed that the polar wind 
mechanism is responsible for the transport of ionospheric H + 
and 0 + , and that only a small fraction of photoelectrons can 
significantly contribute to the ion acceleration. 

In this Letter, we investigate how the ambipolar electric 
field and the fraction of photoelectrons reduce the gravitational 
potential, and therefore, increase the ion scale height and the ion 
density at the top of the atmosphere of HJs. We also calculate 
how the mass-loss rate for H + is modified by this effect. Due 
to the high EUV radiation, the fraction of photoelectron in 
the atmospheres of HJs is expected to be higher than in Earth’s 
case, leading to a much greater increase of electron temperature. 

In Section 2, we calculate the change of the effective gravity 
and the ion scale height due to the ambipolar electric field and 
photoelectrons. We present and discuss the results in Section 3, 
and draw our conclusion in Section 4. 


assuming u a o = u a as a lower limit. Using Equations (2) and (4), 
the electron density at altitude z can now be obtained, assuming 
Zo = 0, r 0 = R p , and r = R p + z: 


n e (z ) = n,(z) - n„(z) = n i0 



? R l \ 

( Rp+z ) 3 / ' 


(5) 


The effective gravity is modified by the ambipolar electric 
field as geff = g — (eE\\/mj), with the ambipolar electric field 
(positive for ions) defined as (Schunk & Nagy 2004): 
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Here T e is the electron temperature, and e is the electric charge. 
At the planetary surface, z = 0 and so we obtain 


2. MODIFICATION OF THE ION SCALE HEIGHT BY THE 
AMBIPOLAR ELECTRIC FIELD 


E,(z = 0) = — 
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In the derivation below, we follow the standard model for the 
polar wind, but we include the effect of the photoelectrons on 
the solution. For a planetary atmosphere consisting of electrons, 
photoelectrons, and ions, charge neutrality requires that 


which yields 


geff = g + 


kT e 
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( 8 ) 


n e o + n a0 = n i0 , (1) 

where n e o, n a o, and n,o are the electron, photoelectron, and 
ion number densities at some reference altitude, r () . From 
Equation (1), we can define the fraction of the photoelectrons, fi, 
as n a o = ( 6 / 1 , 0 , and the fraction of electrons as n e a = (1 — / 6 )/i,o. 

Our goal here is to calculate how the effective gravity at the 
surface is modified when taking into account the photoelectrons 
and the ambipolar electric field, and investigate how this 
modified gravity affects the ion scale height, 77,. We will 
compare 77, with the unchanged scale height Hq which contains 
the surface gravity g but not the ambipolar electric field. 

We begin by assuming a hydrostatic ion density profile: 


n i(z) — riioe {z 4o)/Hi , (2) 

with the ion scale height, Hj = (kf /m,g e ff), where k is the 
Boltzmann constant, 7) is the ion temperature, m, is the ion 
mass, and g e ff is the effective gravity. Without the effects we 
study here, g e ff = g- Conservation of the photoelectron mass 
along a magnetic flux tube requires that 

tla O^GO^O — W a U a A, (3) 

with u a o and u a being the photoelectrons velocities, and 
Aq and A being the magnetic flux tube cross sections at 
the reference altitude and at some altitude, respectively. This 
equation implicitly neglects any scattering of the photoelectrons. 
In a magnetic dipole geometry, the magnetic flux conservation 
requires that AqBq = AB, with B = C/r 3 being the dipole 
field magnitude as a function of radius (C is a constant), and 
Bq = C/tq is the field magnitude at the reference altitude, ro. 
Therefore, Aq/ A = r^/r 3 , and we have 



or 


geff = 


3 fkT e 

mi R p (\ - y6) 


(1 - P)Tj 

T e + ( \ -m 


(9) 


In Equation (9), g is modified by the ion and electron temper- 
atures, and by the fraction of photoelectrons. For the case of 
ft = 0 and 7) = T e , the well-known reduction of the effective 
gravity of the ions by half is obtained (Gombosi 2004). 

As shown by previous models (Tam et al. 1995; Khazanov 
et al. 1997; Tam et al. 1998; Glocer et al. 2012), the electron 
temperature is highly affected by even a very small fraction of 
photoelectrons. In our model here, we assume that 7) = 1000 K. 
Despite the higher ion temperature expected in HJs, the effect 
studied here is driven by the difference between T e and 7), so 
that it should scale with the increase in 7). We scale the electron 
temperature with the percentage of photoelectrons and 7) using 
two different models. One is based on the electron temperature 
distribution at the top of Earth’s atmosphere from Khazanov 
et al. (1997): 

T e (P) = Ti*2 6+logli , (10) 


with 1000 K < T e < 16,000 K for 10- 6 (10“ 4 %) < f < 
10~ 2 (1%), and a more modest function with 1000 K < T e < 
10,000 K: 

T e (f) = f * 1.8 6+log ^. (11) 


With the above models for T e . the modified gravity and scale 
height can be calculated as a function of the fraction of 
photoelectrons. 


3. RESULTS AND DISCUSSION 
3.1. Results 

Figure 1 shows the effective gravity as a function of the 
photoelectron percentage for the two models for T e . For (3 — 0, 
we obtain g e ff/g = 0.5. In Figure 2, we show the electron 
temperature and the ratio of modified to non-modified scale 
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Figure 1. Effective gravity as a function of the photoelectron percentage 
(0.0001%-1%) for T e model 1 (solid line) and model 2 (dashed line). 


height as a function of the fraction of photoelectrons, assuming 
Wi = m p , the proton mass. Here, we show the solution only for 
photoelectron percentages of 0.0001-1. One can see that if the 
fraction of photoelectrons is even less than 1 %, the scale height 
increases by a factor of 2-15. 

A realistic ion density profile cannot be obtained using the 
simplified calculation we present here. In particular, we cannot 
calculate the density profile of the H + ions, since they are 
expected to attain supersonic speeds. Therefore, it is hard to 
estimate the increase in density at the top of the atmosphere and 
the corresponding increase in mass-loss rate. Nevertheless, we 
can use a hydrostatic profile to estimate the ion density change at 
lower altitudes. In Figure 3, we show the ratio of the hydrostatic 
density profiles using the modified and unmodified scale heights, 
respectively, as a function of the fraction of photoelectrons for 
altitudes of 350 km (~ 1 H f} ) and 1000 km (~3//o). The density 
is increased by a factor of 2-3 at 350 km and by a factor of 
5-15 at 1000 km. At higher latitudes, the hydrostatic solution 
is probably not valid any longer and the ratio in Figure 3 will 
become too large, since the density profile for the unmodified 


scale height goes to zero faster than the one with the modified 
scale height. 

3.2. Discussion 

In HJs, the extremely strong radiation is expected to increase 
the fraction of photoelectrons. Therefore, the electron temper- 
ature should be higher than the ion temperature, despite the 
strong heating at the day side, so that the mechanism proposed 
here should still be significant. The effect should be limited at 
the night side due to the lower ionization rate, and it is not yet 
clear how effective the atmospheric day-night circulation is at 
higher latitudes (where the day-night temperature difference is 
smaller than that at the equator), and at high altitudes (where 
the ion acceleration occurs). 

For a magnetized HJ, the mass loss is expected to take place 
along the magnetic field lines which are open to the stellar 
wind (as demonstrated by Stone & Proga 2009; Trammell 
et al. 2011; Adams 2011), and that is exactly where the polar 
wind process takes place. It has been previously shown that 
the classical polar wind mechanism together with the addition 
of photoelectrons and wave-particle interactions is responsible 
for the transport of H + and 0 + out of Earth’s atmosphere. By 
lowering the potential barrier, these processes effectively lower 
the escape velocity. These processes have also been speculated 
to be important at Jupiter and Saturn (Glocer et al. 2007; Nagy 
et al. 1986); the major ions in the upper atmosphere at these 
planets are H 3 and H + . In HJs, the relative ion abundances are 
not known, but modeling by Garcia Munoz (2007) shows that 
H + , HJ , He + , C + , and various ionized hydrocarbons are possibly 
present. The polar wind process should apply to each of these 
planets. Indeed the derivations presented here reflect the basic 
textbook derivation of the classical polar wind (Gombosi 2004) 
to which we have added the effect of photoelectons. No other 
planet specific parameters are required. Even neglecting the 
effect of photoelectrons, the polar wind process by itself could 
significantly increase the ion scale height. 

The relative composition affects the polar wind process by 
changing the parallel electric field. This is because the parallel 
electric field was found to increase with mass. In the case 
of no photoelectrons, if H 3 was the major ion species (such 
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Figure 2. Electron temperature (top) and the ratio of modified to unmodified ion scale height (bottom) as a function of the photoelectron percentage (0.0001%-1%) 
for T e model 1 (solid line) and model 2 (dashed line). 
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Figure 3. Ratio of non-modified to modified hydrostatic ion densities for T e model 1 (solid line) and model 2 (dashed line) at z = 350 km (top) and at z = 1000 km 
(bottom) as a function of the photoelectron percentage (0.0001%-1%). 


as at Jupiter or Saturn) then the parallel electric field would 
exert an upward force approximately equal to one half the 
gravitational force acting on H 3 . In this case the scale height 
of H 3 would increase. Lighter constituents such as H + would 
actually have a net upward force resulting in an eventual 
supersonic flow. Including photoelectrons increases the electric 
field and intensifies the effect of the polar wind, possibly 
resulting in a net upward force on heavier species. If H + was the 
major ion, the parallel electric field would be reduced, but the 
effect would still be quite significant. 

The simplified model presented here is insufficient to predict 
the detailed change in the ion density profile, but it can predict 
how the scale height changes. We show that this change can 
reach about a factor of 10 at lower altitudes. Therefore, it should 
also increase the mass-loss rate by the same amount assuming 
the same surface area and without changing the ion velocity at 
the top of the atmosphere. The polar wind is expected to further 
accelerate the ions such that the ion speed should increase as 
well, so the factor of 10 increase is a lower limit. 

In order to perform a more detailed calculation of the effect 
of the polar wind on the mass-loss rate of HJs, a more detailed 
model is needed, such as the polar wind model by Glocer et al. 
(2007, 2009, 2012), which is similar to that of Garcia Munoz 
(2007), but includes the effect of the ambipolar electric field and 
photoelectrons. The derivation and discussion contained in this 
Letter, however, demonstrate that the polar wind process plays 
an important role in the mass-loss rate of HJs and should be 
accounted for in models. 

4. CONCLUSIONS 

In this Letter, we perform a simplified calculation of the effect 
of the ambipolar electric field and atmospheric photoelectrons 
on the planetary ion scale height. We show that this effect 
can reduce the effective gravity and therefore, enhance the ion 
acceleration in the region of the planetary atmosphere, where 
magnetic field lines are open. We find that a small fraction of 
photoelectrons (less than 1 % of the total electrons) can increase 
the ion scale height by a factor of 2-15. We calculate the 


hydrostatic density profiles using the modified scale heights and 
find that the planetary mass-loss rate should increase by an order 
of magnitude at a minimum, even neglecting any increase in the 
ion velocity due to this the process. Since the ion acceleration 
should be enhanced by the process, we expect the increase 
in mass-loss rate to be even greater. A more comprehensive 
calculation, however, requires a more detailed modeling effort. 

We thank an unknown referee for her/his review report and 
Jeremy Drake for his help in preparing this manuscript. O.C. is 
supported by SI Grand Challenges grant No. 405 10254HH0022. 
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